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Modification of a Variational Objective Analysis Model for New 
Equations for Pressure Gradient and Vertical Velocity in the 
Lover Troposphere and for Spatial Resolution and Accuracy of 
Satellite Data 

1 . Introduction 

Beginning in late 1982, the NASA supported research to 
extend and improve upon the Achtemeier (1975) numerical 
variational objective analysis model to include observations from 
space based platforms for the diagnosis of cyclone scale weather 
systems. The goal of this research is a variational data 
assimilation method that incorporates as dynamical constraints 
the primitive equations for a moist, convectively unstable 
atmosphere and the radiative transfer equation. Variables to be 
adjusted include the three-dimensional vector wind, height, 
temperature, and moisture from rawinsonde data, and cloud-wind 
vectors, moisture, and radiance from satellite data. This 
presents a formidable mathematical problem. In order to 
facilitate thorough analysis of each of the model components, we 
defined four variational models that divide the problem naturally 
according to increasing complexity. The first of these 
variational models (MODEL I), which is the subject of this 
report, contains the two nonlinear horizontal momentum equations, 
the integrated continuity equation, and the hydrostatic equation. 
MODEL 1 is based upon Sasaki's (1958, 1970) method of variational 


objective analysis. Problems associated with an internally 
consistent finite difference method, a nonlinear hybrid 
terrain-following vertical coordinate, formulations for the 
pressure gradient terms, formulations for the velocity tendency 
terms and the development of a convergent solution sequence are 
addressed with MODEL 1 . 

The theoretical development, coding for numerical 
computation, and preliminary testing of MODEL 1 was completed by 
August, 1985. The results of these tests indicated that before 
the theoretical complexity of the model was increased by the 
inclusion of the energy equation as a fifth dynamic constraint, 
it was advisable to modify MODEL 1, if necessary, a) to improve 
the way the large nonmeteorological contributions to the pressure 
gradient force were reduced, b) to generalize the integrated 
continuity equation, and c) to introduce horizontal variation in 
the precision modulus weights for the observations. The results 
of our work on these three topics are the subject of this report. 

Section 2 summarizes our work with the rederivation and 
implementation of an improved hydrostatic equation and pressure 
gradient force algorithm. This work was especially critical to 
the success of the overall MODEL I development after it was found 
that the original derivation was in error. All current and 
future publications including the final report (Achtemeier, 
etal., 1986) present results from the corrected version of MODEL 
1. Section 3 introduces the generalized integrated continuity 
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equation and Section 4 presents the progress toward formulations 
for the horizontal variation in the precision modulus weights. 

2. Reformulation of the Hydrostatic Equation and the Horizontal 
Pressure Gradient Force 

The motivation for changing the method for removing the 
large nonmeteorological contributions by unlevel terrain in the 
hydrostatic equation and in the horizontal pressure gradient 
terms in the nonlinear horizontal momentum equations came about 
from a detailed evaluation of the MODEL 1 variational objective 
analysis model. He had found an anomously large filling of a 
deep synoptic scale trough that was located over the Great Basin 
and Rocky Mountain highlands areas. Our analysis of the 
contributions of the individual constraints to the final fields 
of geopotential height did not reveal that dynamic balancing or 
filtering inherent in the model was the cause for the filling. 
Preliminary calculations indicated that a modification in a 
"terrain correction temperature" could realize a correction in 
the pressure gradient force which, when translated into the 
spatial distribution of geopotential height, was of the magnitude 
of the filling of the synoptic scale trough. Further 

investigation revealed that the algorithm appeared to be correct 
in differential form, for pressure surfaces but there was doubt 
about whether it was valid in difference form on sigma surfaces. 
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Given the uncertainties in the original formulations, it 
was decided to rederive the hydrostatic equation and the pressure 
gradient force terms of the horizontal momentum equations to 
remove the orographic effects in a more rigorous fashion. This 
posed a major reprogramming effort because 1) the terms changed 
are important terms in the variational adjustment and any major 
modifications in their form would result in equally major 
modifications in the Euler-Lagrange equations, 2) the terms are 
expressed explicitly in the derivations to form a diagnostic 
adjustment equation for the geopotential height through the 
reduction of the number of variables and their modification would 
require rederivation of the adjustment equation and introduce 
additional complicating terms, and 3) as described in the 
paragraph below, these terms determine the form in which the 
initial data must enter the variational model. 

The procedure for acquiring the variational ly adjustable 
part of the thermodynamic variables once they have been 
interpolated into the sigma coordinates first requires that all 
variables be nondimensionalized. Second, we remove a hydrostatic 
component that includes much of the vertical variations of the 
lower coordinate surface due to variable topography. Third, a 
hydrostatic reference atmosphere is removed and finally, the 
residual fields are multiplied by the ratio of the Rossby number 
to the Froud number to bring them into compatibility with the 
nondimensionalized dynamic equations. The variational 
adjustments are carried out on these residual variables. 
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Several formulations were made to achieve the requirements 
of the second step. In order to minimize the reprogramming 
tasks, ve sought formulations that did not change the forms of 
the major terms nor introduce nonconstant coefficients. These 
formulations are summarized below. 

a) Mo removal of topographic effects 

The first approach at correcting the formulations for the 
removal of the nonmeteorological contribution to the geopotential 
heights by the underlying coordinate surface was to not remove 
them. The programs for processing the initial data were modified 
so that there was no removal of a terrain component from the 
initial data. Several additional modifications to MODEL 1 were 
necessary to maintain compatibility. Fig. la shows the heights 
of the lower coordinate surface at 1200 GMT 10 April 1979 with 
the large contribution by unlevel terrain included. The figure 
shows mostly the variation of height due to elevation. The 
Appalachian highlands and the Mississippi Valley are visible in 
the east and the western highlands are visible over Colorado and 
Utah. The smoothed elevations, which extend to about 1600 m, do 
not resolve the higher mountain ranges. Fig. lb shows only the 
meteorological heights as they appear on the 1000 mb surface. A 
low about 100 m deep is located over Colorado. This 
meteorological system is an order of magnitude smaller than the 
terrain height. If the latter is not balanced by the pressure 
field which is part of the second term of the pressure gradient 
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fore* as it appears in terrain following coordinates, 
considerable error can be introduced into the adjusted height 
gradients and, through the solution of the dynamic constraints, 
introduced into the vind field. 

Our calculations show that the height field in Fig. la is 
mostly balanced by the pressure field (not shown) through second 
term of the pressure gradient force. However, the variational 
equations separate these terms and build the large 
nonmeteorological heights into the solution. Upon running MODEL 
I, ve found that the cyclical adjustment sequence did not 
converge to a solution. It is imperative therefore that the 
terrain effects be separated from the height adjustments as was 
originally attempted. 

b) Removal of the Nonmeteorological Hydrostatic Component: 
Method A 


Assume that observations of temperature and height have been 
gridded onto the nonlinear sigma coordinate surfaces and 
nondimensionalized according to classical scale theory. Assume 
also that the pressure has been determined in a manner that is 
consistent with the heights and temperatures. The nondimensional 
hydrostatic equation is 


3<t> 


w 


3a 


+ T 


w 


3 In p, 

1 o" 


w 


( 1 ) 


where the subscript, w, identifies the whole or total 
untransformed variable. We partition the nonmeteorological 
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component through the definitions for height and pressure, viz., 

<k* m fa $T 80,1 P w “ P T + p where the subscript, T, 
identifies the terrain contribution. The hydrostatic equation is 
then 


3<J> _ 

-w + ^ T » + 6 ' 


where 


p 3 In p 


'w 


r w 


3o 


6’ 


3<J>, 


_T T w 3p T 

3a P w 3a * 


( 2 ) 


(3) 


(4) 


Equation (2) is identical in fora to the original version in 
MODEL 1. The coefficient JT is horizontally variant however, and 
this increases the complexity of the MODEL I equations over the 
original derivation. 


If we assign the pressure for the lowest three layers to be 
800, 900, and 1000 mb, respectively, then the remaining partition 
variables may be determined by setting ft - 0 and solving for the 
terrain height with the initial, unadjusted temperature for T^. 
Then we remove the reference atmosphere for each sigma level by 
defining <f> , Y, and T w - f R + Y where the 

subscript, R, refers to the reference atmosphere. Substitution 
into (2) leaves 

IT + H + <t r + Y > (t r + T > + 6 ' - 0 * (5) 

We also require that tL be found through 

3* r 

3a + y r t r 0 * (6) 
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Than the hydrostatic equation for the adjustable meteorological 
residual is 

H + V +e - 0 (7) 

vhere ^ • P + Note that though 0 was set to zero in 

order to obtain the terrain height, it will not be zero once the 
intial temperature has been replaced by the variational 
correction. 


Fig. lc shows the heights on the 1000 mb equivalent 
pressure surface after the removal of the terrain contribution by 
Method A. Although the large height anomaly present in Fig. la 
has been removed, the use of Method A is undesirable for the 
following reasons. First, Method A removed too much height from 
the areas of high elevation, ie., it overestimated the depth of 
the low center over Colorado in comparison with the actual 
heights of the 1000 mb surface (Fig. lb). The mean layer 
temperature used for the calculation of the thickness between 
700-1000 mb was estimated by the mean layer temperature between 
700 mb and the pressure of the lowest sigma surface (810 mb at 
the highest elevation). Use of this relatively cold temperature 
in the hypsometric equation underestimates the thickness of the 
layer which, when subtracted from the heights at 700 mb, should 
underestimate the depth of the low at 1000 mb. Thus Method A may 
have overestimated the depth of the low by as much as 120 meters. 
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Second, scale analysis of ths terms of (7) revealed that the 
variational adjustment of the temperature in the second term of 
(4) was equal in magnitude to the correction of the second term 
of (7). Thus (7), written in the form that can be most easily 
adapted into MODEL 1, does not guarantee convergence of the 
solution sequence. Indeed we found convergence only when the 
temperature in the second term of (4) was held left unadjusted. 


c) Removal of the Nonmeteorological Hydrostatic Component: 
Method B 

As in Method A, Method B begins with the nondimcnsional 
hydrostatic equation valid for the whole or total thermodynamic 
variables, 


3 * w _ 3 In p w 

3a + T w 3a 


0 , 


( 8 ) 


We define the meteorological and terrain partitions of the height 
and pressure as before, <j^“ and P w * P<j + P» 

substitute into ( 8) this time expanding the logarithmic term 
also. The resulting hydrostatic equation is 


3<j> 

* TT„ + 8' - 0, 


where 

, „ 3 In p 
3a 


( 10 ) 


p„ . . T W 3p T 


H W . v W . W * T 

s' --sr + ( T ' 15 ^ T ^ ' 

In this formulation, Y is horizontally invariant and thus 


(ID 
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retains the simplified form of the hydrostatic equation as it 
appeared in the original derivation of MODEL 1. We again define 
p - 800, 900, 1000 mb as equivalent pressure surfaces as a 

distinction between the heights obtained by Method B end the 
heights of the original pressure surfaces as obtained directly 
from the observations. If ^ is averaged over each layer to get 
( p R , the reference atmosphere may be removed subject to the 
condition that 


if <J <J^+ ^ and T^ - +T. Then ( 9) reduces to 

H + ^ + 8 - °, 


( 12 ) 


( 13 ) 


where, upon variational adjustments for height and temperature, 
takes the form 


3 

(— - 1 ) — 
'‘P ; 3a 


i 3P X 

0M*>> + 


( 14 ) 


The advantages of Method B are that Y is horizontally 
invariant and that scale analysis reveals that the adjustment 
terms of (14) are at least an order of magnitude smaller than the 
second term of (13). Thus (13) will not impede the convergence 
of the solution sequence of MODEL 1. The difficulty with Method 
B is that the total and equivalent pressures appear as 
undifferentiated coefficients. Thus the accuracy of (14) and 
hence (13) is critically dependent upon a method to determine a 
representative pressure for each sigma layer. After some 
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experimentation, it was found that an accurate representative 
pressure was given by the average of the arithmetic mean plus 
twice the geometric mean, 

P - 3 [j (P 2 + P x ) + 2 1 . (15) 

Fig. Id shows the height of the 1000 mb equivalent pressure 
surface after the application of Method B. The resemblance of 
all features to the heights of the actual 1000 mb surface (Fig. 
lb) is evident except for the smaller central height of the low 
center over Colorado. The underestimation of this feature was 
expected because of the cold layer thicknesses there as described 
in reference to Method A. Since ve have merely partitioned the 
heights, not neglected height terms, the remaining heights that 
make up the difference in the heights between Fig. Id and Fig. 
lb must be contained in the residual term. Method B was adopted 
for the hydrostatic constraint in MODEL 1. 


Once the terrain and meteorological variables have been 
partitioned through the hydrostatic equation, it is a simple 
matter to calculate the pressure gradient force terms of the 
horizontal momentum equations. The nondimens ional pressure 
gradient term in the x-direction is given by 


to 


3 In p 


PGX 


+ T 


w 


3x ' *w 3x * 

* 

where the superscript, x, implies that the temperature has been 
averaged over the interval of differentiation. Partitioning of 


( 16 ) 


the heights and pressures leaves 


PGX 


3x 


3<t> T _ x 3 In p_ 

1 + T X — 5 • 

3x w 3x 


( 17 ) 
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W* not* that before the partitioning, the height of the lover 
coordinate surface was entirely nonmeteoro logical by the 
definition of the terrain-following vertical coordinate. The 
pressure term contained both the terrain and meteorological 
contributions to the pressure gradient force. The partitioning 
has forced part of the meteorological pressure into the heights 
as ^p / ) x *0 and also into through p T in (11). 

Therefore, the pressure gradient is not given totally by the 
first term of (17) as is evident by the differences between the 
height analyses in Figs, lb and Id, but by the first term plus 
elements of the second and third terms. The partitioning is not 
a coordinate transformation, hence the term equivalent pressure 
surface, nor is it a perfect separation of the terrain from the 
meteorological pressure gradient. It is only a method to reduce 
the large nonmeteorological component between the original 
pressure gradient terms. 

Introduction of the reference atmosphere and application of 
scale analysis modifies (17) as follows, 


pgx - + n 

3x x j 


where 


T* 


3 In p. 


3x 


- + R °/r I; <* T + t r ln <>, 


( 18 ) 


( 19 ) 


The temperature in the first term of (19) is subject to the 
variational adjustment. Since the remaining terms are determined 


from initial variables and held fixed, we may rewrite # x aa 


° + (T X -T° X ) 


3 ln p, 


w 


3x 


( 20 ) 
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to giro tho remaining meteorological part of the horizontal 
presaure gradient force. 

3. Generalization of the Integrated Continuity Equation 

The transformation of the continuity equation from pressure 
coordinates into the nonlinear sigma coordinates involved ■ a 
correction term that is the vertical velocity multiplied by a 
coefficient that is a nonlinear function of sigma. This 
coefficient appears as an integral upon the solution of the 
integrated continuity equation for the vertical velocity. An 
exact solution of the integral was not obtained at the time the 
variational model was being coded. Our attempts to approximate a 
solution by averaging failed because of the highly nonlinear 
relationship between pressure and sigma. Rather than hold up the 
model development with attempts to approximate this term, it was 
decided to set the term equal to one and proceed. The result was 
that packing of the sigma coordinate surfaces over elevated 
terrain was not taken into account in the calculation of vertical 
velocities. Low level divergences were accorded too much weight 
in the integration and subsequent adjustment of the continuity 
equation. 

We now know that an exact solution for the integral can be 
obtained if the integration is taken over the pressure interval 
between sigma layers. We thus proceed to generalize the 
integrated continuity equation. The equation of continuity in 
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nonlinear sign* coordinates is 


3u _3v dc 

3x 3y 3cr 


+ 


F 


where 

F - q 2 » 3 



/3u/K L 3v/K x 
' 3x 3y ' 


( 21 ) 


( 22 ) 


In addition, 

-2 [J - ct (p-p*)] 

^ ’ 

and 

(p-p*) 3 J 3 [2ctJ(p-p*)] 

q 2 " (P s -P*) 4 ]2 * 

where 

a - a*/(p*- P(j ) i 


(23) 


(24) 


(25) 


and 

J - 36 (p-p*) 3 + a(p-p*) (26) 

if 

6 • U-o'&j 2 )] <27) 

r *u 

The definitiona of the various symbols are as follows. 

■ dp g /dt is the vertical velocity in pressure coordinates 
along the lower sigma coordinate surface, R^ ■ 0.1, K is the 
variable part of the nondimensionalized Lambert conformal map 
scale factor, p^ ia the surface pressure, p u is the 
pressure at the top of the analysis domain, and p* ■ 700 mb is 
the pressure at a reference sigma level, p*, defined for the 
nonlinear sigma coordinate system. Refer to Achtemeier, etal., 
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19#) for a more complete description of the dynamical 
constraints and the sigma coordinate system. 

When solved for the vertical velocity, (21) takes the 
following form, 


5 ■ 5 o Q o - 0 H Q do , 


(28) 


where 


h - in + Iz + f 

H 3x 3y > 


and 

Q - C q l *' 


(29) 


(30) 


We seek the solution of the integral term Q. Upon conversion to 
an integral over pressure, the exponent of (30) becomes 


P 

* n 


q L da 


1 3cr 
q l 3p dp • 


(31) 


From the definition of the sigma coordinate, 


a » f5(p-p*) 3 + a(p-p u ) j 

the derivative, dff*/dp, is found to be equal to 
- 3B(p-p*) 2 + a , 

Combining (33) with q^ in (23) and integrating over pressure 
gives 

’ P 1 q |2L d p - - In [3B(p-p*) 2 + 

P 

and therefore, (30) becomes 

3B(p 1 - P*) 2+ct 

Q - 2 ' 

3B(p-p*) +ot 


(32) 


(33) 


(34) 


(35) 
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Because the term, Q» i* a variable coefficient that 
multiplies the horizontal divergence* it causes the introduction 
of additional terns in the Euler- Lagrange equations for the 
variational adjustment model. These terms increase the 
complexity of the diagnostic equation for the adjustment velocity 
potential. Therefore* the Euler-Lagrange equations which are 
affected by the generalization of (28) are rederived. All new 
terms are introduced into the equations in a manner consistent 
with the vertically staggard grid of the variational model. Then 
the variational model is rerun snd a series of tests conducted to 
verify the coding of the formulation and determine how the 
correction term impacts upon the vertical velocity. 

The third constraint* M^, expressed in complete form is 


M 


3 



(u +v ) e 
x y 



q x da 


da + (a-a 

O 



da 


) 


+ f 1 >1 ’2 + *1 ^ '“ A 1 




- \ (u'V y + v’V*) ] e 


«1 da * 


(36) 


Performing the variation only on the product of the Lagrangian 
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multiplier with (36) gives the following terns in the adjustment 
equations for u, v, and (T. 


6 u [X 3 M 3 3 



0 Mo) * - \ « 3 f 1 0 *** 

O °o 


do) 

x 



da) 

) 


(37) 


5v [X 3 M 3 ] 



p Q da) y - ^ (X 3 J^ 1 Q K* 7 dff) y 


" h (X 3 


C q k * w j 

o 


(38) 


6o [X 3 M 3 J * X 3 # 


(39) 


Upon inspection of the terms of these equations, it becomes 
apparent that neither X 3 nor the integral of Q is a function 
of sigma. Therefore we can combine both into a new Lagrangian 
multiplier, 



o 


The variations on u, v, and < 7 * (eqs. (37), (3 8 ), and (39)) 
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simplify to 


5u [X 3 M 3 ] - - X 3 X ~ ®i <^ 3 kX7 ) x " *i (X 3^c ) , 

y 

5v [X 3 M 3 ] * “ x J y “ *1 (^^ icy > y ~ (X 3 K y^ ) 


( 41 ) 

( 42 ) 


6or [X 3 M 3 ] 



( 43 ) 


These modifications are inserted into the full 
Euler-Lagrange equations that contain reference to ^ 3 . The 
Euler-Lagrange equations for u and v are 


V - X 3* + X 2 + F 1 


( 44 ) 


F 1 * - n i u ° + E 1 C> ' X 2 + R o 



- lm* 5 'x F (u-c x ) IC l j[ - [m (v-c y ) x J y 


— x-xya 


- ^ (oV ) } - ^ I(X^) x + x 3 T? 1 


L W 


( 45 ) 


v- 



+ F„ 


0 

> 


( 46 ) 


_ o _ — x 

- IIjV - c 


\ + \ ♦ ^ 7 % 
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( 47 ) 


— y--xya 

- [m X 2 (u-c x ) 7 ) x - [m^X* (v- c y ) 7 ] y - R q (a X ) a > 


- R-L t(X^K Xy ) y + X^Kj ] . 

The Euler-Lagrange equation for 0* i* 



( 48 ) 


TV * O 

-It 2 a 




— X 

u a 



( 49 ) 


In addition, the integrated continuity equation ia 



Q(u x +v y ) da + a + F 7 


0 


( 50 ) 


q o a° + 




- ^ (u*K^ + ] Q da 


( 51 ) 
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Variable* are eliminated to produce a diagnostic equation in 
> Dividing (44) and (46) by 77J and reformulating as 
components of the divergence gives 


U x A 3 x J x + [ n x (X 2 + F l )3 x 


(52) 


V 

X 



3yV + 



F 2»y 


0 . 


(53) 


Then (52) and (53) combine into the divergence and after 
integration over some interval dtf, become 


1 * 
Q(u +v ) da + (A 
x y 3x 


p* da) x 

Ja o n i 


+ (A 


3y J 


1 d°>. 


a n, 

o 1 


/ + F l )J x - \ ( - X l + F 2»y> « d11 ■ 0 ' 

0 


(54) 


Nov (f is eliminated from (49) and (50) so that 
* 

t0 -\ *3 F 

Q(u +v ) do - - — r— 3 + F 

* y n «> °i F 7 

1 Q do u 2 


0 . 


A diagnostic equation in A 3 is obtained upon elimination of the 

integrated divergence through combining (54) and (55): 

* 


(A 


3x 


°1 0 * r ° 

; \ 

O 0 


f°l £ 
1 *1 


da). 


Oo 

* n 


+ F 


10 


Q da 


where 


(55) 


(56) 


10 


1 <>2 + V'x + >1, ( - X l + F 2»y } q ** - ^ + F 7 • 


(57) 
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We also note that several terns of (56) obey an identity (Eq.(48) 
in Achtemeier et al . , 1986). Therefore (56) transforms into the 
two-dimensional second-order elliptic partial differential 
equation with non-constant coefficients given by 


t<S 


1 da X 
U 1 


* 

3xx 


+ 


fC l o y * 

1 e-) do x, + 

o n i 3 > 1 

o 


0l (S-) 

o n i x 
o 


do X 


* 

3x 


f 1 \ 


) do X, 
y 3y 




-ft + F 


10 


Q do 


(58) 


In the event that jT t is horizontally invariant, (58) further 
simplifies to 


ro 

'a 


1 

0 


(S-) do v 2 x 
U 1 


* 

3 



+ F 


10 


do 


0 


(59) 


Model 1 was rerun vith Q calculated from the above 
development and compared with Q"1.0 for all levels. Fig. 2 
shows the vertical velocities at level 3 (700 mb). There is 

found an approximate 50 percent reduction in the vertical 
velocity for the analysis where variable layer thickness is a 
factor, the reduction having been made over the higher elevation 
areas west of the Great Plains. For example, a greater than 4 cm 
sec ^ center of rising motion for the Q*1.0 analysis (Fig. 
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for 


2a) has baan reduced to slightly greater than 2 cm sec 
the analysis with variable Q (Fig 2b). Elsewhere over lower 
elevations the impact of the layer thickness is much smaller. 

4. Impact of Horizontal Variation of Precision Moduli 

As a simplifying factor in the derivation of the MODEL 1 
variational analysis, we made the precision moduli that weight 
the observations functions of height only. This is a reasonable 
assumption as long as the data are relatively evenly distributed 
over the analysis domain and the data quality are horizontally 
invariant. The TIROS-N data violates the assumption on both 

counts. The T1R0S-H data we are using for the 10 April 197 9 
model verification case studies contains large gaps which are of 
regional scale. The method for objective gridding of the data 
fills the gaps by interpolation from surrounding data. The 
resulting temperature fields are only approximate in these areas 
and can be significantly in error if the data gaps coincide with 
temperature maxima or minima. Our analysis of 12 GMT 10 April 
197 9 was such a case. 

Model 1 was designed for horizontally variant weights. Ho 
rederivations were necessary for this part of the study. We 
began sensitivity studies with the velocity adjustment potential 
equation number (58). The precision modulus weights were allowed 
to vary over five orders of magnitude giving adjustments that 
ranged from 100 percent restoration of the horizontal velocity 
divergence to 100 percent restoration of the initial vertical 
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velocity. The weights were distributed over a number of 
geometrical patterns in order to determine the conditions which 
violate the stability criteria for (58). The results to date are 
too preliminary to report upon at this time however, early 
indications are that, in some of the cases we tested, the use of 
the five order of magnitude variation in the weights, even though 
satisfying stability criteria, can produce greatly exaggerated 
horizontal velocity fields. 
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Figure 1. Heights of the lower coordinate surface for MODEL 1 at 
1200 GMT 10 April 1979, a)the unlevel terrain included, b) the 
meteorological heights as they appear on the 1000 mb surface, c) 
the heights on the 1000 mb equivalent pressure surface after the 
removal of the terrain contribution by Method A, and d) the 
height of the 1000 mb equivalent pressure surface after the 
application of Method B. 
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Figure 2. The vertical velocities at level 3 (700 tab), 1200 GMT 
10 April 197 9 for a) for the Q-1.0 analysis and b) for the 
analysis with variable Q. 
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